Effects of interface spin-orbit coupling on tunneling between normal metal and chiral 

p-wave superconductor 
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We study the tunneling conductance of a clean normal metal/chiral p-wave superconductor junc- 
tion using the extended Blonder-Tinkham-Klapwijk formalism. It is shown that the spin-orbit 
coupling of the Rashba type that is present near the interface causes the subgap conductance peaks 
associated with the Andreev surface bound states to shift to a nonzero bias. We also investigate the 
effect of the Fermi wavevector mismatch between the normal metal and the superconductor. 
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I. INTRODUCTION 

Tunneling spectroscopy is one of the most powerful 
probes of the electronic states in superconductors. Quasi- 
particles in anisotropically-paired superconductors can 
experience a variation, e.g. a sign change, of the or- 
der parameter upon reflection from an interface. Then 
the interference between the incident and reflected quasi- 
particles results in the formation of bound states near 
the interface, with the energies inside the bulk energy 
gap, which are known as the Andreev bound states 
(ABS),^'^ The ABS manifest themselves in low-energy 
features, typically peaks, in the tunneling conductance, 
which have been observed experimentally. In particu- 
lar, the zero-bias conductance peaks in high-T^ cuprates 
exhibit strong dependence on the crystallographic orien- 
tation of the interface, consistent with the d-wave pair- 
ing, see Ref. ^3 for a review. The observation of broad 
subgap peaks in the tunneling conductance of Sr2Ru04 
(Ref. Ij) can also be explained in terms of the sur- 
face ABS, which are expected to exist in a chiral p- 
wave superconductor i^i^ Other systems studied recently 
include the interfaces between magnetic or nonmag- 
netic normal metals and noncentrosymmetric or mag- 
netic superconductorsi^"— We would like also to mention 
that the zero-bias conductance peaks can also originate 
from the quasiparticle states localized near strong impu- 
rities or surface inhomogeneities.-^ Those can be distin- 
guished from the ABS by their different response on a 
magnetic field iii 

Due to the breaking of reflection symmetry near the in- 
terface, quasiparticles experience the spin-orbit coupling 
(SOC) of the Rashba type^^ even if both the normal 
and superconducting crystals have inversion symmetry 
in the bulk. The effects of such interface SOC have been 
neglected in the previous studies of the tunneling con- 
ductance. In this paper we focus on the properties of 
a junction between a normal metal and a chiral p-wave 
superconductor. It is known that the ABS in this case 
correspond to Majorana fermions with linear dispersion, 
see Refs. IgIIisI and [ij- There are strong experimen- 
tal indications that the chiral p-wave state is realized in 
Sr2Ru04 (Ref. Il5l ). We neglect disorder and calculate 



the tunneling conductance using the Blonder-Tinkham- 
Klapwijk (BTK) formalismJ^ The effect of the Rashba 
SOC on the conductance can be attributed to a modifi- 
cation of the boundary conditions for the wavefunctions 
at the interface. We also take into account the differ- 
ence between the Fermi wavevectors on the normal and 
superconducting sides. 

The paper is organized as follows: In Sec. II, 
we develop a theoretical model of the normal metal- 
superconductor (N-S) junction with the interface SOC 
and use the BTK approach to calculate the amplitudes 
for various quasiparticle scattering processes. In Sec. 
Ill, the effects of both the interface SOC and the Fermi 
wavevector mismatch (FWM) on the tunneling conduc- 
tance are presented and discussed. Sec. IV contains a 
summary of our results. In Appendix, we analyze the 
ABS spectrum in a half-infinite chiral p-wave supercon- 
ductor with an arbitrary interface potential. Throughout 
the paper we use the units in which h — \. 



II. FORMULATION OF THE MODEL 

We consider a two-dimensional clean N-S junction 
shown in Fig. [TJ The interface is located at a; = 
and is characterized microscopically by a potential bar- 
rier which is not necessarily an even function of x, due to 
the crystal structure difference between the normal and 
superconducting sides. The asymmetric part of the po- 
tential is responsible for the interface SOC of the Rashba 
type. We consider the following model for the interface 
barrier: 



U{x) = [Uo -t- Uin ■ [a- X k)]5{x), 



(1) 



where n = i; is the unit vector along the interface normal, 
C/q and f/i are the strengths of the the spin-independent 
and the Rashba SOC contributions, respectively, a are 
the Pauli matrices, and k = — iV. The band disper- 
sions are assumed to be parabolic, with the same effective 
masses of quasiparticles on both sides (according to Ref. 
IitI . the effect of the mass difference is equivalent to that 
caused by a variation of the interface potential strength) . 
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FIG. 1: (Color online) Schematic illustration of the quasipar- 
ticle reflection and transmission processes at the N-S inter- 
face. 



On the superconducting side of the junction, we assume a 
chiral p-wave pairing state of the form d{k) ex z{kx+iky). 
The Bogoliubov-de Gennes (BdG) equations for the 
four components of the quasiparticle wavefunction are 
decoupled into two independent pairs of two-component 
equations as follows: 



■H^*(r)=£;*(r), 
where a — ± for different spin orientations, 



(2) 



Ho 



Zm 



2m 



At(fe,r) T^ + Ef; - U„{x) 



(3) 

and Ua{x) — (Uq — aUiky)5(x). The Fermi energies in 
the normal and superconducting regions can be differ- 
ent due to different carrier densities, with Epi = EpN or 
Eps- The ratio of the corresponding wavevectors is char- 
acterized the dimen sionless FW M parameter as follows: 
Ao = kps/kpN = \/Efs/Epm- 

The off-diagonal elements of Eq. (j3]) contain the gap 
function A. In the spirit of the BTK approach, we neglect 
self-consistency and assume the gap magnitude to be 
equal to Aq on the superconducting side, and zero on the 
normal side. Then we have A(fc,r) — {Ao/2kps){{kx + 
iky),9{x)}, where ^(a;) is the step function. The anti- 
commutator on the right-hand side is required since the 
order parameter varies in space, see Ref. |lj. Below we 
use a simpler expression: A = (Ao/fcFs)(fca; + iky), at 
a; > 0, neglecting the (5- function term in the off-diagonal 
elements, which gives a small correction to the boundary 
conditions. 



Suppose an electron is injected from the normal metal 
with the excitation energy E > and spin cr, at an angle 
6 from the interface normal. The momentum parallel to 
the interface is conserved: 



kpN sin 6 = kps sin Og 



(4) 



The incident electron is reflected back either as an 
electron (normal reflection) or as a hole (Andreev 
reflection) ji^ In the superconductor, the wavefunctions 
of the transmitted quasiparticles have both electron and 
hole components. 

Solution of Eq. (g]) has the form *(r) = e'''yy^{x), 
where 



*w(a;) = ( ' ^ 



ikfN COS 9 X 



ik^N COS 9 X 



+ ba 



-ikpN COS 9 X 



(5) 



on the normal side, and 
^s{x) = c 



ve 



r^ikps COS 6s X 



-ve 
u 



-iOs 



-ikps cos 6s X 



(6) 



on the superconducting side. Here a„ and b„ are the 
amplitudes of the Andreev and normal reflection, respec- 
tively, and Co- and d„ are the transmission amplitudes. 
The electron and hole components of the wavefunctions 
in the superconducting region are given by 



1 / Q. 

T2r^r 



1 / Q. 



(7) 



where Q. = ^J E^ ~ l^l. 

All the reflection and transmission amplitudes in Eqs. 
([5|) and (jni) can be found from the boundary conditions 
that follow from Eq. ([T]): 



*S(0+) = *Ar(0-), 

*s(0+)-n(0~) 

= 2m(L/o - cr[/ifc_FArsin6l)*jv(0"). 



(8) 



In particular, for the reflection amplitudes we obtain: 
aa{E,e) 



AXuj^e-'" 



baiE,0) 



(f + A2 + Z^)uj+ + 2Acj_ ' 

[(1 - iZ^f - \^]uj+ 
(1 + A2 + Z^)uj+ + 2Acj_ ' 



(9) 



where 



A = Ao 

u 

V 



cost 



COS 9 ' 

UJQ = -, UJ+ 
V 

Zo 



cos 6* 



aZi 



sin 6* 



-2i0 



-2i6 



COS 6*' 



Zo 



2mUo 
kpN 



Zi = 2mUi. 



The dimensionless parameters Zq and Zi characterize the 
strengths of the potential and SO scattering, respectively. 
From the conservation of the probability current it fol- 
lows that 



+ A(|c„ 



K 



1. 



We note that the charge transmission at subgap energies 
is enhanced when the normal reflection is minimized, i.e. 



J 



at 6cr = 0, which happens if uj+ = 0. Then Eq. ^ 
yields a subgap resonance with the energy E = Aq sin 9s , 
ds > 0, corresponding to a branch of excitations localized 
near the interface. This result agrees with that of Refs. 
lallSl. andH. see also Appendix. 

By using the BTK formalismji^ we obtain for the di- 
mensionless angle-resolved differential tunneling conduc- 
tance: 



G^{E,9) = l + |a,(^,0)|2-|6,(^,0)|2 

^ 4A{[(1 -I- A)^ + Z'M + 4Ac.g - [(1 - A)^ + Zg]} 

[(1 + A)2 + Z2]2c^4 + 2COS20[(1 + A)2 + Z2][(l - A)2 + Z2]cj2 + [(^ _ ^)2 + ^2p 



(10) 



r 



One can see that, while the conductance depends on the 
incident spin orientation: G+ {E,0) ^ G- {E,9), the time 
reversal invariance is respected: Ga{E,9) = G-a{E,—9). 
The experimentally measurable tunneling conductance 
G{E) is obtained after the angular integration and the 
summation over the incident spin orientations as follows: 

1 ,^ /•'r/2 
G{E)^-^Y d9cos9Ga{E,9). (11) 

Gn „ J-t,/2 



the purely potential and the SOC contributions to the 
interface scattering. In the top panel, we show G{E) in 
the absence of the interface SOC {Zi = 0). The mid- 
dle and bottom panels demonstrate the effects of varying 
the interface SOC at Zq = (high transparency barrier) 
and Zq = I (low transparency barrier). Note that, if 
Zq = Zi = 0, then G{E) = 2 for all subgap energies, 
-E < Ao, due to the Andreev reflection occuring with the 
probability one. 



Here Gat is the conductance for a normal metal/normal 
metal junction with the interface potential given by Eq. 
©: 



with 



Gat^V/ d9cos9GNa{E,9), 



Gn(j{E,( 



ZI 



(12) 



(13) 



We note that the difference between the Fermi wavevec- 
tors imposes a constraint on the effective range of an- 
gles contributing to the integrals in Eqs. pT|) and 
(|12p : if kps < kpjsf, then there is no transmission for 
I sin 6*1 > Aq. 



The broad peak at the subgap energies, which is most 
pronounced for a low-transparency interface, is associ- 
ated with the surface ABS that exist in chiral p-wave 
superconductors^ (we recall that there are no subgap sur- 
face bound states in s-wave superconductors, and the 
tunneling conductance at £' < |A| is strongly suppressed 
for large Zq, Ref. [la ). When we include the interface 
SOC, the conductance peak is shifted to a nonzero bias, 
see the bottom panel of Fig. 2. The origin of this effect 
can be attributed to a nonmonotonic angular dependence 
of the transmission coefficients at Zi ^ 0, which is ev- 
ident from Eq. ([T5|) . Note that the spectrum of the 
surface ABS is not affected by the interface SOC, see 
Appendix. 



III. RESULTS 

The tunneling conductance of the N-S junction at zero 
temperature can be plotted as a function of the dimen- 
sionless excitation energy E/Aq. We will study the ef- 
fects of the interface SOC and the difference between the 
Fermi wavevectors on the tunneling spectra. 

Let us first consider the case in which there is no FWM, 
i.e. kpN = kps and Aq = 1. Fig. 2 shows the tunnel- 
ing conductance G{E) for different relative strengths of 



Next, we discuss the effect of the FWM on the tun- 
neling conductance, in two cases: (i) Aq — 1.2, (ii) 
Ao = 0.7. In each case, we consider both high- and low- 
transparency interfaces, Zq — and Zq = 1. Figs. 3 and 
4 show the variation of G{E) for several values of Zi . In 
the case (ii) (and for kps < kpN, in general), the con- 
ductance is notably suppressed by the FWM, because 
the modes with \sm.9\ > Aq experience total reflection 
and, therefore, do not contribute to the conductance. In 
contrast, there is no discernible consistent effect of the 
FWM in the case (i). 
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FIG. 3: The dimensionless tunneling conductance G{E) ver- 
sus E/Ao for the FWM parameter Ao = 1.2 and difTerent 
values of the interface SOC. In the top panel Zq — (no po- 
tential barrier), in the bottom panel Zq = 1 (strong potential 
barrier) . 
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FIG. 2: The dimensionless tunneling conductance G{E) ver- 
sus E/Ao for Ao = 1 (no FWM) and different values of Zq and 
Zi: Zi = (top panel), Zq — (middle panel), and Zq = 1 
(bottom panel). 



IV. SUMMARY 

In summary, we have applied the extended BTK for- 
malism to investigate the tunneling conductance of a 
junction between a normal metal and a chiral p-wave 
superconductor. We focused on the effects of the inter- 
face SOC of the Rashba type, both with and without 
the Fermi surface mismatch between the two sides of the 
junction. The structure of the subgap peaks in the tun- 
neling conductance has been shown to strongly depend on 
the interface SOC. In particular, for a low-transparency 
interface, the maximum of the conductance associated 
with the surface ABS is shifted away from the zero bias. 
We also considered the case of different Fermi wavevec- 
tors in the normal and superconducting regions. When 
kps < kpN, the tunneling conductance is strongly sup- 
pressed. 
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Appendix A: Surface ABS Avith spin-orbit coupling 

Let us consider a superconductor with a flat surface, 
occupying the a; > half-space. We assume a hard-wall 
confining potential V{x), which is infinite at a; < and 
varies within a thin surface layer, whose thickness a is 
of the order of several lattice spacings. This potential 
results in the SOC of the form (1/4to^c^)F'(x)(p x &)x, 
which is also restricted to the vicinity of the surface (re- 
call that ft = 1 in our units). Assuming a 2D geometry 
with the band dispersion ^(fe) (which includes the chem- 
ical potential), the quasiparticle spectrum can be found 
from the following 4x4 BdG Hamiltonian: 



n = 



A 



At 



(Al) 



where e — ^(— iV) + V{x) — 
single-particle Hamiltonian. 



ii/4m^c^)V'{x)a3\'y is the 
In a triplet pairing state, 
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FIG. 4: The dimensionless tunneling conductance G(E) ver- 
sus E/Ao for the FWM parameter Ao — 0.7. Other parameter 
values are the same as in Fig. 3. 



the order parameter is given by A = (ia'a2)d{k,r). To 
make analytical progress, we neglect self-consistency and 
assume that the order parameter is uniform at all x > a. 
The system is translationally invariant along the surface 
and the wavefunctions have the form ^(r) — e''"'!'^^(x). 
For a chiralp-wave state of the form d oc z{k,j;+iky), it 
is easy to show that Eq. (lAip can be written as a direct 
sum of two 2x2 Hamiltonians given by 



Ho 



_ fC + UAx) 



At 



A 

-i-u^{x) 



(j{l/Am^c^)V'{x)ky, 



(A2) 



and the 



a = ±, UAx) = V{x) 

gap function is given by A = {Ao/kp){kx + iky), with 

kp = kps being the Fermi wavevector. 

In the region x > a, Ua{x) — 0, and the spectra of 
%(, can be analyzed in the semiclassical, or Andreev, 
approximation^ We represent wavefunctions in the form 






"(/"(a;), where kx satisfies the equation 



(A3) 



at fixed ky. Each Fermi-surface wavevector A; — {kx,ky) 
defines a semiclassical trajectory, along which the quasi- 
particle state is described by a coherent superposi- 
tion of the electron and hole amplitudes: ip{x) = 



\%pe{x),il^h{x)Y' , satisfying the Andreev equation: 



-iVF,x{k)^x ^(fe) 

A'*(fc) ivFAk)^x 



i:{x) = Eijjix), (A4) 



where vpik) = d£,/dk is the quasiparticle velocity on the 
Fermi surface, £' > is the energy of excitations, and 



A(fc) = Ao^^^^^. 



(A5) 



Depending on the direction of propagation, the semiclas- 
sical trajectories are classified as either incident (vp^^ < 
0) or reflected {vp^x > 0). 

One can seek solution of Eq. ()A4|) in the form of a 
plane wave: ij{x) ~ e"''^. Focusing on the quasiparticle 
states which are bound to the surface, but cannot exist in 
the bulk, we expect that E < Aq. For the wavefunction 
we then obtain: 



^k{x) = 



= ±f I 



where 

a{k) 



V2 V "C^) 
A*(fc) 



E + iVt{k) sigT\vp^x{k) 



^ — K.{k)x 



K{k) 



(A6) 



n{k) 

\vF,x{k)\' 



and n{k) = ^|A(fc)|2-£;2 = ^A^ - E^. 

The Andreev approximation is not valid near the sur- 
face, at < X < a, where the rapidly varying potential 
Ua{x) is nonzero. The surface scattering will result in 
the effective boundary conditions at a; = a, which ex- 
press the Andreev wavefunctions corresponding to the 
reflected trajectories in terms of those corresponding to 
the incident trajectories.— 

Depending on the band structure and the surface ori- 
entation, Eq. (jA3|) might have several solutions. In the 
case of a parabolic band, there is only one incident and 
one reflected trajectory at each ky (except ky = ±kps, 
where the Andreev approximation is not applicable): 
kin = {-K,ky), kout = {kx,ky), with k^ > 0. The 
effective boundary condition has the following form: 



V'feo„t(a) = Sipk,„{a), 



(A7) 



where S is the surface scattering matrix, which is an 
electron- hole scalar— determined by the surface potential 
U„{x). Inserting here the wavefunctions (|A6p . we arrive 
at an equation for the bound state energy: 



E + in{k,n) _ A{kin) 



(A8) 



E-in{kout) A{kout) 
For the chiral p-wave state (jA5l) , the solution is given by 



E{ky) = AoT^ {ky > 0). 



(A9) 



Remarkably, this expression does not contain any micro- 
scopic details (in particular, it has the same form for "H^ 



and H-), which is consistent with a topological nature 
of the surface ABSi^ We come to the conclusion that 
the surface bound states of the Hamiltonian (lAll) are de- 



scribed by two {a = ±) degenerate branches of fermionic 
excitations with linear dispersion, given by Eq. (jA9p . 
regardless of the SOC strength. 
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